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Abstract 

A parabolic stochastic PDE is studied analytically and numerically, when a bi- 
furcation parameter is slowly increased through its critical value. The aim is to 
understand the effect of noise on delayed bifurcations in systems with spatial de- 
grees of freedom. Realisations of the nonautonomous stochastic PDE remain near 
the unstable configuration for a long time after the bifurcation parameter passes 
through its critical value, then jump to a new configuration. The effect of the non- 
linearity is to freeze in the spatial structure formed from the noise near the critical 
value. 
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1. Introduction 

Many physical systems undergo a transition from a spatially uniform state to one 
of lower symmetry. Such systems are commonly modelled by a simple differential 
equation which has a bifurcation parameter with a critical value at which there is 
an exchange of stability between steady states [1]. The model can be improved 
by adding noise [2], which can be understood as allowing for small rapidly varying 
effects neglected in the formulation of the model. Noise is also necessary to pro- 
vide the initial symmetry-breaking which permits the system to choose one of the 
available lower-symmetry states. 

If the bifurcation parameter is not a function of time, the effect of noise is to 
make a bifurcation ill-defined in an 0{e) region near the critical value, where e is the 
size of the noise [3] . A more dramatic effect is found if the bifurcation parameter is 
a function of time, corresponding to an experimental situation in which the critical 
value is not known in advance, so the parameter is slowly changed until a qualitative 
change is seen in the behaviour of the system [4]. 

In the case of the pitchfork bifurcation, described by the normal form 



it is normally said that, if g > 0, y will be found either at y = ^/g or y = — y^. 
However, on solving the stochastic differential equation 



with g = fit, starting with y = at ^ = < 0, it is found that y remains near 
y = until well after ^ = (Figure 1). (The Wiener process - standard Brownian 
motion - is denoted by lower case w.) If /i is small, a typical trajectory consists of 
a long sojourn near y = 0, a sudden jump away from y = and then relaxation 
towards = g. The value of g at which the jump occurs is a random variable whose 
probability distribution can be found by solving the linearised version of (2) [5]; the 
mean value is approximately A/2//[hie[ and the width of the probability distribution 
is proportional to fj,. 

The phenomenon of delayed bifurcation just described and its sensitivity to 
noise has been studied theoretically and experimentally in the case of non-autonomous 
stochastic ordinary differential equations [6-9]; the corresponding phenomenon, de- 
layed transition, for PDEs is the subject of this article. A dynamic transition is 
described by the following non-autonomous parabolic stochastic partial differential 
equation: 



y^gy- y^ 



(1) 



dy = {gy — y^)dt + edw 




dY = {g{t)Y -Y^ + AY)dt + edW, 



(3) 
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Figure 1. Dynamic pitchfork bifurcation with noise. The solid line is one 

trajectory of the ordinary SDE (2), with t = lO"^'^, /x = 0.01, and initial condition 
y — 0.1 at g = —0.8. Also shown, as dotted lines, are the loci of stable fixed 
points of the corresponding autonomous system (1). The trajectory lingers near 
y = for a long time after g passes through 0; the value of g at which the jump 
towards y ~ 'i^^/g occurs is a random variable with mean approximately equal to 
-^2/i| lne\. 

where g = litis slowly increased through 0, W is the Brownian sheet, A = Ym!=i 

i 

is the Laplacian in R"* (the derivatives exist only in the sense of generalised func- 
tions) and jJL and e are constants such that 



(4) 



In section 2 this initial value problem is solved using periodic boundary con- 
ditions in one space dimension, in which case y is a stochastic process with values 
in the space of real- valued continuous functions on [0, L], and has continuous mean 
(Yt{x)) and correlation function (Yt{x)Yt{x')) [10, 11]. In section 3 the correspond- 
ing equation in space dimension, m, greater than one is discussed. 

The present work was motivated by the behaviour of some nonlinear systems 
of ordinary and partial differential equations whose dynamics are controlled by 
noise [12, 13]. The controlling influence of noise arises because trajectories spend 
long times near a slow invariant manifold, and the dynamics near this manifold are 
similar to those of a dynamic bifurcation [5] . The most important parameter in the 
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Figure 2. Dynamic transition in one space dimension. One 

numerically-generated realisation of the stochastic PDE (3) Y-i{x), x G [0, L], 
is displayed at four times as the bifurcation parameter, g = f^t, is slowly increased, 
passing through at t = 0. (Note the different vertical scales.) It remains close 
to the zero configuration until well after g = 0. Nonlinear terms become important 
when g ~ -y/ 2fi\ In e|; their effect is to freeze in the spatial structure formed, 
near g = 0, from the noise. The data shown here were generated using a finite 
difference method with 500 points. (L = 300, /x = 0.01, e = 10" ^^.j 



description of these systems is n\ lne|, where ^ is the timescale for the dynamics 
near the slow manifold and e is the magnitude of the noise. 
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2. Model dynamic transition in one space dimension 

Four configurations from a numerical reafisation of (3) with m = 1 are shown in 
Figure 2. According to static stabihty analysis the zero configuration is unstable 
for > 0, but Y typically remains much closer to than to = ^ until well after 
g = when g is time-dependent. The first part of a dynamic transition, the loss 
of stability of the zero configuration, is thus controlled by the noise and the sweep 
rate, not by the nonlinearity. 

If Y obeys (3), the linearised version of the stochastic differential equation for 
its kth Fourier mode, it^^-* = ^ Yt{x')e'-'^^^'^/ ^dx' , is 

du^^^ = {g{t) - {^fk^)dt + edw^^\ (5) 
ij 

with each w^^^ an independent Wiener process. Thus, for as long as nonlinear terms 
are unimportant, each Fourier mode evolves independently, satisfying an equation 
of the dynamic bifurcation type. 

As the in ODE case, the new configuration (Figure 2(d)) appears abruptly at 
a value of g which is a random variable with mean approximately lne|. The 

effect of the cubic nonlinearity is to arrest the exponential growth of (Yt{x)Yt{x)) 
when it becomes 0{1) and to freeze in the spatial structure formed during the slow 
sweep. Establishing when this happens makes it possible to answer another question 
which is not considered in the static version: what is the typical size of the spatial 
domains formed in the transition? 

To examine the dynamics of the loss of stability of the zero configuration, 
consider the solution of the linearised version of (3), 

dY = {g{t)Y + AY)dt + edW, (6) 
with initial data Yt{x) = f{x) on a; e [0, L] ai t = to < when g = go: 

Y{x,t)= f G{t,to,x,x')f{x')dx' + e f [ G{t,to,x,x')dWt{x'). (7) 

Jo Jto J[0,L] 

The first integral is an ordinary integral; the second a space-time Ito integral. The 
fundamental solution G{t, to, x, x') is given by 



I 1/2 /2\ ^ ^ {x — x —jL) 

G(t,t',x,x')= , e^^(^ ) y e . (8) 

^ ^ V47r(t - 1') ^ ^ ' 



j = -oo 



The solution (7) prescribes the mean and correlation function of y as a function 
of time. The mean value of Y is the first term, an ordinary integral which, for large 
t — to, is given by 

(rt(x)> = e^'^(*'-*o)y (9) 
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where j/ is a smoothed version of the initial data f(x). In the rest of this article, it 
is assumed that y < 0{1) and 2/i|lne| < Qq. This ensures that (Y^) <^ (Xt)i ^® 
Y is taken as a mean-zero Gaussian process. It is also assumed that |to| > so 
that quasi-static equilibrium is reached before 5^ = 0, and L ^ \/t — to so that, at 
any one x, only one term in (8) is important. 

The correlation function is a mean square quantity which can be evaluated 
using the 'delta function' property of the Brownian sheet as an integrator: 



{Yt{x)Yt{x')) = e^ f [ G^{t,s,x,x')dsdx' 

Jo Jto 

= ds , =Q 8(t-«) . 

Jto V87r(t - s) 

For times t such that —t > l/y/ju the nonautonomous correlation function (10) 
differs only by C(/x) from the result obtained for fixed g < if to ^ oo [14]: 

{Yt{x)Yt{x')) = ^^e-l"--'IVl^. (11) 

In the non-autonomous case, the correlation function remains finite as t passes 
through 0, and for large positive times {t > 1/-^), it is well approximated by: 

{Y,{x)Y,{x')) (12) 

Note the Gaussian form of the spatial correlation and the existence of a character- 
istic length y/si at time t. This formula is valid until g ~ •\/2/i| lne|, when the cubic 
nonlinearity becomes important. The correlation length at this time becomes the 
characteristic size of the spatial domains formed (Figure 2(d)). 
3. Higher space dimensions 

Proceeding as in section 2, the expression obtained for the correlation function at 
time t is 

{Ytix)Ytix')) = / ds—- -^e~^^, (13) 

Jto (87r(t-s))2 

where ||a; — a;'|| is the Euclidean norm on R"^. The major difference from the 
case m = 1 is that the correlation function docs not approach a finite limit as 
||a; — a;'|| 0. This is a general feature of second order stochastic PDEs in more 
than one space dimension [15, 16]. On a finite grid, equation (6), with g < fixed, 
is of interest because it generates a stochastic process with non-zero correlation 
length in both space and time [17]. The singularity in the correlation function at 
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X = x' manifests itself in the normalisation, which increases without a finite limit as 
the grid spacing is decreased. In the nonautonomous version the same divergence 
is seen (logarithmic when m = 2), but for > the singularity is only noticeable 
on very small length scales. The nature of these divergences, in both discrete and 
continuous versions, is conveniently studied in Fourier space, where the linearised 
equation separates into uncoupled ordinary SDEs, each of which can be solved 
exactly. 

The finite difference method for a parabolic SPDE consists of replacing the 
infinite dimensional system (3) by A^'" ordinary SDEs on a grid of equally-spaced 
points in [0, Lf^. For the SPDE (3), the SDE at each point is 

dY{x) = {gY{x) -Y^{x))dt + pAY{x)dt + p'^edw (14) 
where the discrete Laplacian A is defined by 

AY{x) = ^(^') - 2mF(x) (15) 

x' 

and the sum is over the 2m nearest neighbours of x. Each of the N'^ ordinary 
SDEs has an associated Wiener process, independent of all the others, and so can 
be numerically solved in the normal way [18]. Note the scaling of the magnitude 
of the noise added at each grid point with p. The smallest length scale resolved is 
roughly the distance between nearest neighbours, As in numerical solution of 
deterministic parabolic PDEs, there is a maximum timestep, proportional to ^, if 
the finite difference method is to be numerically stable [19]. 

Solved on a finite grid in this way, the behaviour of realisations of the non- 
autonomous version SPDE (3), with m > 1 and g = iit as in section 2, is qualita- 
tively similar to that in one space dimension [20]. That is, a Gaussian form of the 
spatial correlation function emerges from the slow sweep through = 0, 




and metastable domains of positive and negative Y are formed when the nonlinearity 
becomes important, at ~ •\/2//| lne|. 
4. Conclusion 

In systems with spatial degrees of freedom, the phenomenon of delayed bifurca- 
tion is accompanied by a noise-controlled formation of spatial structure. Here the 
stochastic PDE (3) was studied, whose realisations consist long sojourns near the 
zero configuration during which the shape of the spatial correlation function be- 
comes Gaussian, followed by an abrupt change to the nonlinear regime, with the 
formation of domains of positive and negative Y . This spatial structure, formed 
near the critical value of the bifurcation parameter g when noise dominates, is frozen 
in at ~ -^2^1 lne|, when the characteristic length is -\/8t- 
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